NPS ARCHIVE 



1968 

LLOYD, R. 


— 


• . 


ANALOG SIMULATION OFAUTOMATIC GLIDE 
SLOPE CONTROL USING LIFT SPOILERS 
AS DIRECT LIFT CONTROL 

by 

Robert Collins Lloyd 
and 

James Kenneth Swift 



library ^ 

NAVAL POSTGRADUATE SCHOOL 
MONT ■ ' ^ 93940 



DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 

MONTEREY. CA9394MiJl^ 



UNITED STATES 

NAVAL POSTGRADUATE SCHOOL 




THESIS 



ANALOG SIMULATION OF AUTOMATIC GLIDE SLOPE CONTROL 
USING WING LIFT SPOILERS AS DIRECT LIFT CONTROL 

by 

Robert Collins Lloyd 
and 

James Kenneth Swift 







f 



4 





ANALOG SIMULATION OF AUTOMATIC GLIDE SLOPE CONTROL 



USING WING LIFT SPOILERS AS DIRECT LIFT CONTROL 



by 



Robert Collins Lloyd 
Lieutenant, United States Navy 
B.S., Texas Technological College, 1960 



and 



James Kenneth Swift 
Lieutenant, United states Navy 
B.S., University of Detroit, 1960 



Submitted in partial fulfillment of the 
requirements for the degree of 

AERONAUTICAL ENGINEER 

from the 

NAVAL POSTGRADUATE SCHOOL 
June 1968 



mAeciitv^ 

U-OHO, A. 




ABSTRACT 



The use of wing lift spoilers as a means of changing lift 
without changing angle of attack was studied for use in the 
landing approach task. The vehicle used was the F-8 type fighter. 
Automatic glide slope controllers were proposed using an elevator 
glide slope coupler in conjunction with an automatic power 
compensator for comparison with an automatic direct lift control 
system. The system gains were optimized for gust disturbances 
and initial offsets from glide slope. An analog computer simulation 
program including a manual control phase was used to determine 
arbitrary measures of effectiveness of the proposed systems. 
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CHAPTER I 



INTRODUCTION 

Since low aspect ratio, swept wing airplanes exhibit a 
relative insensitivity to change in lift with angle of attack 
and a sharp increase in drag with increasing angle of attack, 
a means of varying lift at constant angle of attack is desirable. 
Such a method is called Direct Lift Control (DLC) because the 
airplane is given a near instantaneous vertical acceleration, 
whereas changing lift by elevator control has an inherent time 
lag followed by an overshoot due to the airplane's moment of 
inertia about the pitch axis. 

Various forms of DLC are: 

1. Lift jets 

2. Rotary wings 

3. Fast acting wing flaps 

4. Wing lift spoilers 

Forms 1 and 2 are not applicable here. Form 3, fast acting wing 
flaps, is treated in References 1 through 5. 

Form 4, wing lift spoilers, was first used as DLC on assault 
gliders in WW II (6). More recently, Bray and Drinkwater of NASA 
Ames Research Center have flight tested wing lift spoilers as 
DLC in a large subsonic jet transport; their work is awaiting 
publication. It is expected that Boeing Airplane Company will 
employ spoiler DLC on the B2707 SST. The implementation of wing 
spoilers as DLC is discussed in Chapter II. 
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At first glance, spoilers seem a negative approach to the 
problem of close glide slope control; however, one must note 
that an overall increase in the level of lift is not the intent 
for which DLC is employed in this context. What is sought here 
is rapid and precise control of the level of lift in order to 
maintain an instrument approach glide slope. Lift control is 
accomplished by trimming the airplane with spoilers extended to 
some base position; then retraction or further extension from 
this datum produces the required change in lift. An increase of 
stall speed is expected when the airplane is trimmed with the 
spoilers extended. The resulting increase in approach speed 
may be offset by use of an automatic power compensator which 
allows lower approach speeds. Also, more precise control of 
speed enables the approach to be made at a shallower glide slope 
angle which reduces the rate of descent. Objectionable airframe 
buffeting caused by extended spoilers may be reduced by venting 
the spoilers. Of great importance is the favorable drag change 
when a lift increase is commanded. This is not the case when 
wing flaps are used as DLC and is certainly not true when elevator 
control is used. This decrease in drag with increase in lift 
shows its full importance in the discussion of automatic systems 
in Chapter V. 

In anticipation of Category III operations (zero ceiling- 
zero visibility landings), very precise glide slope sensors and 
transmitters are being designed. The Navy has recently developed 
an Automatic Carrier Landing System (ACLS) (7) which can bring 
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airplanes aboard ship in zero-zero weather. Obviously tight glide 
slope control is of paramount importance under these circumstances. 

Automatic glide slope control, discussed in Chapter III, is needed 
to provide this capability. 

The problem of automatic glide slope control is twofold; glide 
slope must be maintained and airspeed held constant. Automatic 
control of airspeed through the use of an Automatic Power Compensator 
(APC) is discussed in References 6, 8, 9, and 10. Automatic glide 
slope control using DLC with and without APC, Pitch Control System 
(PCS) without APC, and an Elevator Glide Slope Coupler (EGSC) with 
APC are proposed in Chapter III. 

The linearized, longitudinal equations of motion for small 
perturbations were used in conjunction with concepts from elementary 
feedback control theory. Numerical values of the stability deriva- 
tives for the Chance Vought F-8 from Reference 9 are used in the 
equations of motion and are listed in Table I. Loop gains for 
the systems were approximated by the root locus method using a 
digital computer program. The gains thus obtained were optimized 
for various disturbances using the analog components of a Comcor 
Incorporated Ci 5000 Hybrid Computer. 

In Chapter IV a man replaced the automatic systems in the 
control loop. Potentiometers were affixed to an aircraft type 
stick and throttle to provide control input signals. Real time 
display of deviation from trim airspeed and programmed glide slope 
was accomplished through a single channel 12 inch oscilloscope. A 
time sharing program was used to present two independent traces 
s imultaneous ly . 
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The effectiveness of the various systems was determined 
from the simulated response of the airplane to horizontal and 
vertical gusts and to an initial glide slope deviation of 10 
feet. Discussion of the results appears in Chapter V. A 
summary of conclusions and recommendations is found in Chapter VI. 
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CHAPTER II 



DLC IMPLEMENTATION 



General 

Spoilers were used for DLC in this study because of their 
near instantaneous response and favorable drag characteristics. 
A flow separation device such as a spoiler is impossible to 
analyze using potential flow theory. There are, however, 
empirical means presented in the literature for the prediction 
of rolling effectiveness of spoilers in high speed flight (11, 
12, 13). When these methods are extended beyond their limits 
to the high angles of attack used by swept wing aircraft in 
the landing approach, the results vary greatly. The predicted 
spoiler lift effectiveness from various sources is as follows: 
Reference 11, Q = 1.4 

i ^ o 

Ref erence 12 , C =1.6 

Ref erence 13 , Q =2.1 



Reference 14, 







= 2.4 



Spoiler Analysis 

Since analysis was next to impossible and empirical methods 
yielded inconsistent results, wind tunnel data were relied upon. 

A search of available data revealed a case which closely approximated 
the design in question. Reference (14) includes a remarkably similar 
wing planform to that of the F-8 aircraft. See Figures 1 and 2. The 
model was full size and was tested with a fuselage in place. 
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In the approach configuration, the angle of attack, > 
range is 11 to 16 degrees for the F-8 (15). Hence the angle 
of attack range for DLC operation was established. 

Pilots use vertical accelerations of order + O.lg to 
maintain a given glide slope (15). Therefore, the system was 
called upon to provide ^(^^of + 0.1 in the DLC range. 

Reference 16 states that spoilers have little or no lag 
in operation when located at the 0.7 c position. Since ailerons 
at or near the wing tip are used for lateral control, the spoiler 
must not extend past the 0.6 b/2 location. The following con- 
figuration provides the required and meets the above position 

constraints . 

X,= 0.7c N;-0.Zb/7 No- 0^6 b/2 



For this configuration, (3,^ vsQ{^curves were plotted from data 
in Reference 14 and appear in Figure 3 for various spoiler deflections. 
Spoiler deflection, 3 > measured in units of non-dimensional 
spoiler projection normal to the wing surface, Vls/c • sense 

of S considered positive- for spoiler retraction since this 

action causes an increase in lift. Spoiler lift effectiveness , , 
is a positive quantity in that an increase in lift is caused by 
positive spoiler deflection. This is not the case when the 
system is used. Conversion to the X* Z system is shown in Table I. 

Changes in constant for various 3 plotted in 

Figure 4. The curves of 0( Figure 4 do not pass through 

the origin. The effect was attributed to the re-attachment of the 
boundary layer when spoiler extensions are small. 
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curve 



The slopes of ZiCi.vs 0( in the I 2 .5 OC ^ 1 4- . <o ^ 
range yielded - 2 . 4 - . Applying to the ClVsO( 

for the F-8 (15), Figure 5, gave the required in the DLC 
range. 

Drag curves (14) are shown in Figure 6. It was noted that 
virtually all of the drag in the DLC range is due to(X, and there 
is surprisingly little change due to spoiler deflection. Since 
the drag change due to spoiler deflection is so small, an expanded 
scale plot of Q^vs (X.for the DLC range is given in Figure 7. 

Spoiler drag effectiveness, 0^^ ’ computed from the/^(^^vs g 

curve in Figure 8. For 0(=/4° > • Results of applying 

Coj, to the F-8 Co vs curve (15) are shown in Figure 9 where the 
steep increase in drag with (X is readily apparent. 

The pitching moment curves are shown in Figure 10. Negative 
static stability of the NACA model in the \7° range 

was due to the absence of a horizontal tail. When the effect of 
the tail is added, the system is statically stable at all 
shown in the F-8 (Xcurve. 

Equations of Motion 

The effect of DLC on the dynamics of the airplane is shown 
in the longitudinal equations of motion. The controls fixed case 
was applicable here because of the assumption of power operated, 
irreversible controls . 

In anticipation of real time analog flight simulation where 
a human pilot would be used, the equations of motion were used in 
dimensional form. All time derivatives were taken with respect 
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to real time as opposed to non-dimensional time which is 
common practice in stability and control work. 

The dimensional equations of m.otion for the longitudinal, 
controls fixed case were taken from Reference 17 and altered 
to include DLC terms. The equations are: 



where the axis system and angle convention is shown in Figure 
11. The quantities on the right hand sides of Equations 1 are 
first order terms of a Taylor*s series expansion. The spoiler 



equations to represent an input from DLC. A suitable elevator- 
DLC interconnect was assumed to compensate for trim changes due 



Cl'- + UXwCX - gCos 0 9 ^ 7^ + Xt i - 



IDOL'S- g 

0= MaU^Ufiwa*UMwa ^ 



( 1 ) 




added to the basic 



to spoiler deflection; therefore, was not introduced into 

the control equations. 
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CHAPTER III 



AUTOMATIC GLIDE SLOPE CONTROL 
I. SYSTEM DESCRIPTION 



APC 

In order to provide speed stability for the EGSC and to 
increase the speed stability of the DLC system, an APC was 
incorporated in this study. The APC investigated in Reference 
9 Was used. This system incorporates an automatic throttle 
controlled by feedbacks of UL and (^ . In anticipation of 
restrictions on the number of operational amplifiers available 
on the Ci 5000 analog computer, the system was modified slightly. 

All time delays were deleted except for the engine acceleration 
time lag. 

Automatic DLC 

An automatic DLC system was considered without APC in order 
to investigate the ability of the system to maintain a given 
glide slope without the artificial speed stability supplied by 
the APC system. This step was prompted by the favorable drag 
characteristics mentioned in Chapter I. The automatic DLC 
controller used incorporates both position and rate feedback. 

Automatic DLC and APC 

The automatic DLC was coupled with the APC in order to determine 
if system performance could be improved by increasing the speed 
stability of the automatic DLC. 
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EGSC and APC 



The conventional method of controlling attitude, and thus 
glide path, is with the elevator. The EGSC was selected to 
provide a basis for the evaluation of the DLC systems mentioned 
above. As shown in Chapter V, the automatic control of glide 
slope with an EGSC is impossible without some form of artificial 
speed stability. 



II. SYSTEM ANALYSIS 



APC 

The block diagram for the modified APC system is shown in 
Figure 12. With the exception of system time lags, the controller 
is the same as the one described in Reference 9. The speed of the 
aircraft is sampled and compared with the desired approach speed. 

If an error exists, a variation in thrust is commanded to eliminate 
the error. A parallel control loop samples variations in (X and 
commands a thrust variation in a similar manner. These thrust 
variations are summed and fed into the airplane aerodynamics. The 
blocks indicate individual transfer functions which will be derived 
later. Standard block diagram algebra (18) was used to obtain the 
APC system transfer function. The result was: 

LL - Ll^ T 

Uc' , ^XiJLX _ i X oc ■ (2) 

' To Tear 

Automatic DLC 

The block diagram for the automatic DLC system is shown in 
Figure 13. The system consists of an outer loop which incorporates 
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position control with position feedback. The position controller 
itself consists of a proportional control and an integration term 
which is supplied to eliminate steady state error. The integral 
term may be thought of as memory since its effect is to make the 
actions of the controller depend upon the history of the error. 

The restoring force is proportional to the product of the average 
value of the error and time. If a small error continues to exist, 
the restoring force continues to increase with time. 

The inner loop of the automatic DLC system incorporates rate, 
or derivative, feedback. Rate feedback has the effect of giving 
the controller the ability to anticipate errors and thus increase 
the effectiveness of the controller. Thus the output signal of 
the spoiler controller depends upon both position error and the 
rate at which position is changing. 

An acceleration type feedback was considered for the inner 
loop but was discarded in favor of the rate feedback system because 
of the roughness encountered by higher derivative controls. 

As in the case of the APC, standard biock diagram algebra 
was used to obtain the overall system transfer function. The 
result was : 

h - hg S ^ V\ (3) 

he I + 

he ^ ^ V\ b h 

EGSC 

The block diagram for the EGSC is shown in Figure 14. The 
controller is identical to that used for the automatic DLC system 
except that inputs to the airplane aerodynamics are elevator 
deflections instead of spoiler deflections. 
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The closed loop transfer function for the EGSC is, from an 



analysis similar to the above: 



1 XU h 

±L- he n y ^ 

he" 1 ii Jl _ 

he>J ^ h ^ h 



(4) 



Transfer Functions 

Basically a transfer function is the ratio of the Laplace 
transforms of the output of a system to the input. The overall 
transfer functions of the systems, Equations 2, 3, and 4, are 
made up of the individual component transfer functions. 

In the block diagrams of the various glide slope control 
systems and the APC, the airframe can be thought of as a plant 
which produces LI > CX > 0 inputs of and ^ or . The 

airframe transfer functions of interest in this study were: 



T 



, and where is the glide slope angle 



perturbation. 

In order to derive the airframe transfer functions. Equations 
1 were recast into matrix form and use was made of Laplace trans- 
form notation. Equations 1 are then: 





“UXvj 9COS0 




~u 




'X, 




X,' 




u 


( ^ " Zw ') ^ q >5inQ 

^ \3 




a 


- 


2x 

u 


U 


■ 2^1 
u ' 


S 




-udiw-id^) . 


1 


.0. 




n, 


Hi 


n,. 
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where the forcing functions are grouped on the right hand side. 



The individual functions were obtained by 


Cramer's rule. 


For 


example, ^ , 


the 


airframe * s 


response in 


airspeed to 


a change 


in thrust at 


constant control deflection 


is : 










X, 


-ux. 


gCos0 










u 


- Zw 








a _ 




Mt 








(5) 


T 






-UXw 


gCos© 








_Za 

u 




.^43^ 










-hv. 








where the denominator is the 


characteristic equation 


(CE) 


of the 


airframe and 


is common to all airframe transfer functions 


. The 


others are: 








9 Cos 0 
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~ Z.UV. 
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-A 4qSin0 

^ 3 u 




(6) 
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( 8 ) 



4 -X. -ux„ X8 




e 

s 



- Mw M% 



Ct 



But 






(9) 



hence 




( 10 ) 



The reader is reminded that the above relations are but four 
of the component transfer functions in Equations 2, 3, and 4, 
Although an analog computer was ultimately used to determine gains 
for the systems, an estimate of the magnitude of the gains was 
necessary so that an iteration process could be used. Accordingly, 
certain well known assumptions were used in the simplification of 
Equations 5, 6, 7, and 8, 

The purpose of the u loop in the APC is to control long period 
oscillations of airspeed. The well known assumption that the 
longitudinal motion of the airplane can be separated into long 
period and short period oscillations was employed here. Long 
period oscillations in u were assumed to occur at constant Cy . 
Equations 1 then become: 




(lb) 
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thus 



T "0 



(5a) 



The (X loop in the APC is provided so that drag due to angle 
of attack perturbation is compensated by thrust. In accordance 
with the assumption above, variation in (X occurs in the short 
period oscillation where u is assumed constant. For the short 
period case Equations 1 become: 






cX 

o 




(Ic) 



thus 

_a ^ 

The automatic DLC controls short period deviations from 
the glide slope, hence, the constant airspeed assumption is used 
for the DLC transfer function. Equations 1 for the DLC under the 
short period assumption become 



A-Zvj 



-A 




(Id) 



thus 
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Rate of deviation from glide slope is given by Reference 3 
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hence 
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From Equation 12 

h-^u^at 

then 
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The remaining component transfer functions are given in Reference 
9 and presented here for clarity. The controller transfer functions 
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relate thrust commands to deviations in C( or u from desired 
values. The u loop controller transfer function is 



Jc 

Ue 






term is the integral term discussed above. Ko. 



where the 

and are the gains which must be determined. TheCXloop 



controller transfer function is 

We 



K. 



where another undetermined gain. The engine is approximated 



by a first order time lag 



T 






ic 

where spin up lag for the engine and is assumed here to 

be 1.15 seconds (9). Other lags such as servo lags and airspeed 
and angle of attack sensing lags were neglected for the reason 
mentioned earlier. 

The remaining component transfer functions in the DLC over- 

C u 

all transfer function, Equation 3 and Figure 13, are -y- , 

0 n 



and 



h 



The position controller, referred to as the glide slope 



deviation controller, has the following transfer function 

8 






Ki 



where the term is provided so that steady state errors in 
glide slope are eliminated. |<^y^and are gains to be determined. 
The deviation rate controller transfer function is 

Sh _ . 

TT ■ 
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where Kv^ is another gain to be determined later • The rate of 
deviation from glide slope is integrated to produce glide slope 
error. The error in glide slope is fed back to the glide slope 
deviation controller which commands a corrective spoiler 
deflection. This spoiler deflection command is summed with the 
command from the rate controller and then fed into the airplane 
aerodynamics . 

The component transfer functions for the EGSC system are 
analagous to those for the automatic DLC system. 

Estimation of Gains 

In each of the foregoing systems there are three unknown gains 
to be determined. Due to the relative ease with which system 
parameters can be varied on the analog computer, its use was highly 
desirable in the selection of loop gains. As stated earlier a 
reasonable estimation of the gains was needed as a basis for gain 
optimization. Since the inner loop on each system involved only 
one unknown gain, the inner loops were analyzed independently and 
gains were selected for a damping ratio of = 0.8. This procedure 
is customary in control engineering for the analysis of multiloop 
systems. The selected gains were inserted into the individual 
transfer functions, and the inner loops were incorporated into the 
overall transfer function for each system. The integration gain, 

Kl , was neglected in the overall transfer function analysis since 
its presence rather than its value is important in the steady state 
operation of the system. The outer loops were analyzed and gains 
were selected as close as possible to y = 0.5. The reason for 
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selecting a higher damping ratio for the inner loop was that 
addition of the outer loop tends to decrease inner loop 
damping. A more complete analysis than the one shown below 
can be found in the Appendix. 

The transfer function for the inner loop of the APC system 



is 



Tc 



T 



T ■ 

Tc CXe T 

Introducing Equations 6a, 16, and 17 along with the numerical values 
of the stability derivatives from Table I the transfer function 
becomes 



^%O.S06i 1 1281 

T' I.Z8\ ■>4,55*10“' K, 



in which the engine time lag has been neglected to allow a simple 
second order analysis. The characteristic equation for the inner 
loop is 

^NO.aOfei) *\.281 *4!)5*IO'‘’K,-0 

from which 

K=<-- -148,500 Ib/rad 

Using this value for Ko< , the inner loop transfer function becomes 

X - I 28 \ 

T ' ,0.806a <0.65 

With the inner loop transfer function determined, the APC 
transfer function was found to be 



ULc" 







35 



Inserting the transfer functions and stability derivative values 
as above, the APC transfer function is 



y=- , jlO.OOI45_Ka./* O.OOII7 K>.^_->_OjOOI86KuJ _ 

^4(d8&6»0,dbl45Kj4’ ’(0734 O,00117K0A’'(a07’Q00iK3j 

A root locus plot was made for the above system using a digital 
computer program. The value for Ky found from Figure 15 was 

Ku.= 4Z I b 5ec/ft 

The APC transfer function then becomes 

a . 0.Q^ ;lSQ.0494tn.07ft .i 

ur 4‘^0.9Z6i’*O7834 tO.l466fO.OZ3 

The gains thus selected were used as initial values in the gain 

optimization procedure, 

•1 

In a similar manner, the inner loop transfer function for the 
automatic DLC system was found to be 

Vi_ 85 4( 0.5 6 4S0. 1 54^ t 0.40 

8k ' f34(O.8Ofo3843K07, '(1281-31.4 K43 -%Kv> 

Root locus analysis of this transfer function, Figure 16, yielded 

K(.= -0.003 sec/U 

from which 

-h = £ 34(0.56^^0.1.54/' ,0 41 ) 

Sh 1.06y tl.5s^ + 0.2 8 3 

Analysis of the outer loop of the automatic DLC system then yielded 

h . Z 54K»(Q.564tO. 1.54 4 r O. 4 I 

nO 4* t (.06 , (\.38t84.3Kv,)44'(0.288t5l.4RO,w96k'K 
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From a root locus plot of the above system, Figure 17, 

Kv, = 0.005 

The overall transfer function for the automatic DLC system is 

h - 0.4 09/^% 0.157^ ^ 0.48 

ho' +0,48 

Since the EGSC system is useless without automatic speed 
control, and since the analysis of the EGSC system coupled with 
the APC system would involve a seventh order transfer function, 
it was decided to attempt gain optimization for this system 
without estimated values of the gains. The iteration procedure 
involved is discussed in the next section. Gain Optimization. 

In order to demonstrate the lack of feasibility of conven- 
tional elevator without APC as an automatic glide slope control- 
ler, a simple Pitch Control System (PCS) was devised. Figure 18. 




FIGURE 18 

BLOCK DIAGRAM-PCS SYSTEM 



The inner loop transfer function was found to be 

_Q. . -(2.25i . 0. B67 ^ 

^ “ ^“■.(0,80b-2,25KJji(l.Z8T-Q8G7Ko^ 

from which, for ^ = 0.8, 

Kg = -O.G67 sec 
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Inserting this v^alue, the overall transfer function was found 
to be 



Gain Optimization 

Within the framework of the feasible gains de.terrained above, 
it is possible to find the best values of the gains tu give optimum 
response to various disturbances. To this end a systematic 
iteration process was employed through the use of the Ci 5000 
analog computer and a fast response machine plotter. 

The disturbances used were: 1) initial glide slope error of 

10 feet low (fly up command), 2) 5 knot tail gust, and 3) 5 knot 
up gust. The reasons for choosing these disturbances were the 
difficulties encountered when flying into the turbulent wake of 
the island structure and flight deck of an aircraft carrier (7). 
Figure 19 shows the geometry of the mirror approach and the relative 
location of the airplane when wake turbulence is encountered. The 
10 foot error corresponds to a full low deflection on the mirror 
landing system at 800 feet from touch down. The gusts and glide 
slope deviations are of necessity small so that the small pertur- 
bation assumption is not violated. Relatively small disturbances, 
however, can produce divergent oscillations if the systems are 
unstable. It was decided that the glide slope systems would be 



e -A2.2.5K^i lOMlM 

e<,' ^ (l,85S-Z.25Ke')^-0.867R9 



From the root locus plot of this transfer function 

K© ’ - 1.3 rad/sec 
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more sensitive to an h command disturbance, while the APC would 
be susceptible to a tail gust. The plan was then to optimize 
the systems for the critical disturbance while maintaining stable 
response to the other disturbances. 

The airplane and the systems were patched into the Ci 5000 
according to Figures 20 through 24. Individual potentiometer 
settings are listed in Table II. The scaling equations are shown 
in Table III. Gains for the APC, DLC, and EGSC were initially set 
as determined from the root locus analysis. 

In the above analysis of the APC, the integration gain, Ki, 
was assiamed zero. For optimization purposes, a small value of Ki 
on the order of .01, was used as an initial setting. A feasible 
value of 1<^ from the analysis was set; then various values of the 

Ol 

u loop gain,|^^^, were tried until a stable system with an acceptable 
rise time of less than, say, 4 seconds was determined. This inter- 
mediate value of Ku was then held constant. Then, was varied 

until an intermediate value of determined. Finally, the 

product was varied until steady error was minimized. With 

*12, Ku was varied about its intermediate value using 
overshoot, rise time, and damping as criteria for choosing a of 
400 as optimum. A similar procedure using overshoot as a standard 
optimized )<^^at -10,000. The gains thus obtained are quite different 
from those predicted by the analysis. This was expected since the 
loops were assumed independent which they obviously were not. How- 
ever, a starting point was all that was desired from the analysis. 
Figures 25, 26, and 27 each show three values of the gains used in 
r>ptimizing the APC system. 
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The procedure for optimizing the gains for the DLC system 

was similar to the above. Values of Kh> KkKl 
.005, -.0073, and iQ'Sespec tivaly , were det.rmintd. Figures ^8, 

29, and 30 show the effect on the system of varying the gains fre^ 
the optimum. Of particular interest is Figur- 28 which confirtr.s 
what was said earlier, that increasing outer loop gain, Kh> decreases 
the damping of the whole, system. Figure 29 illustrates a divergent 
response when the rate loop gain, Kh> decreased in an effort 
to obtain a shorter rise time. As expected, the integration gain 
had little effect on the damping of the system. 

Gain optimization for the EGSC APC combination was performed 
directly on the analog computer. It was found that a gross esti- 
mation of the gains was all that was necessary if a few more 
iterations were used. Thus the somewhat lengthy root locus analysis 
was avoided in this case. Optimum gains for > Kh Kh 

found to be \. 15 \\ 6 ^ > 4.4)^10^ . 

Effects of varying loop gains are depicted on Figures 31, 32, and 
33. An interesting effect is noted on Figure 32 where an increase 
in K,:, produces separate motion superimposed on the basic oscillation. 
The source of this motion was not determined, but it was noted that 
the trace became more erratic as was increased. This effect 
established an upper boundary on the inner loop gain. When over- 
shoot is limited to, say, 20 per cent of the initial displacement, 
the rise times for the EGSC APC combination were twice as long as 
rise times for the LLC system. This was the first indication of 
the superiority of the DLC system over the EGSC APC combination. 
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CHAPTER IV 



ANALOG SIMULATION PROGRAM 



General 

The Analog Simulation was accomplished in three phases* 

The first phase involved the analog solution of Equations 1 to 
establish the responses of the basic airplane without control 
as a standard. Automatic control was introduced in phase two. 

In the third phase a man was placed in the control loop. 

The Comcor Inc. Ci 5000 analog computer was used in this 
study. A photograph of the control console is provided in 
Figure 34. The computer has 52 operational amplifiers installed 
with space provided for 84 more. Of these 52 amplifiers, 20 
may be used as integrators. The machine incorporates 48 servo 
set potentiometers along with 32 manual pots. The servo setting 
feature allows rapid changing of system parameters. Time scaling 
for the speeding up or slowing down of the solution is readily 
accomplished on the logic patch board. The analog computer 
combined with a SDS 930 digital computer comprise the NPGS Hybrid 
Computer. 

Basic Airplane 

The linearized longitudinal equations of motion for small 
perturbations were used in the Analog Simulation. Numerical values 
of the stability derivatives for the F-8 were substituted in Equations 
1. The equations were then scaled for the 100 volt Ci 5000 according 
to the scaling equations and scale factors in Table IV. After 
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substitution, scaling and rearrangement, Equations 1 appear as: 



a =-OJOia-O.OU(;>et-0.!l ie-0.252T-O.00Gy -0.0055 
a -Q544a-Q43a- 5 - 0.01 935 G-0 082?y-0.202§ 
§ =0.031 a-0.398C(-0.0l5ft-0. 1 16 6 -0.7857^-0.091 T 

where indicates scaled variables , 

The. analog circuit diagram for the basic airplane is seen 
in Figure 20. All such diagrams used herein employ standard 
symbology; the reader is referred to any good text on ana.log 
computer techniques such as References 17, 18, or 19. Although 
the diagrams are not as concise as they might be, it was felt 
that clarity was more important than style. Potentiometer 
settings are found in Table II. 

In phase one no control was provided; hence the basic 
airplane's inherent dynamic stability was relied upon to close 
the loop. The basic airplane was perturbed from the trimmed 
state by the following gust disturbances: 

1) Tail gust of 5 kt magnitude 

2) Up gust of 5 kt 

3) Combined 5 kt tail and 5 kt up gust 

As stated earlier these disturbances are reprc-sentative of 
conditions astern an aircraft carrier making 35 knots in still air. 
Due to separated flow from the angled deck and the island structure 
there is turbulence and a defect in velocity extending about 800 
to 1000 feet aft of the ship (7, 15). The reader is referred again 
to Figure 19. 
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The basic airplane was also subjected to step control inputs 
of stabilator and DLC spoiler in order to compare their effects. 

Automatic control systems introduced in phase two were: 

1) APC, 2) DLC, 3) DLC plus APC, 4) EGSC plus APC, and 5) PCS. 
Circuit diagrams for each component system are provided in Figures 
21 through 24. Scale factors are found in Table III, and 
potentiometer settings are in Table II. 

The APC is used to augment the speed stability of the basic 
airplane and cannot be considered by itself as a means for glide 
slope control. Therefore, the APC system was subjected to the 
same gust disturbances as was the basic airplane. In this case 
the task for the APC was to maintain trimmed airspeed in the 
influence of a gust. 

The glide slope controllers were required to maintain glide 
slope in the influence of gusts. Additional required tasks were 
that of returning the airplane to the glide slope from an initial 
offset of 10 feet low and the combination of initial 10 feet low 
offset and 5 knots slow. 

The PCS system was checked only for its response to a step 
pitch command. Its inclusion was for the purpose of showing the 
infeasibility of using automatic elevator without an APC for glide 
slope control. 

Phase three of the simulation used a man in the control loop 
assisted in some cases by the APC. The pilot was provided with a 
stick and throttle (Figure 35). Visual cues were provided by 
traces on a 12 inch oscilloscope pictured in Figure 34. Photographs 
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of the analog and logic patching used for the manned simulation 
are shown in Figure 36, 

Potentiometers were affixed to the stick and throttle through 
gear trains to assure adequate amplification of the small control 
movements anticipated. Ten turn potentiometers required an 
elaborate transmission to effectively use the resolution they 
provided; hence, the use of these was discarded in favor of one 
turn pots with a simple gear train. The gear trains and pots 
were located at the base of the stick and inside the throttle 
assembly . 

Artificial feel was provided by centering springs on the 
stick. The throttle assembly included a friction control to 
suit the pilot *s preference. 

At first, visual display was tried using the multiple channel 
oscilloscope which is part of the analog computer accessories. 

It was found that by increasing the time scale, oscillations in 
pitch were displayed as a vertically translating horizontal line 
analogous to the attitude gyro in an airplane. It was not possible 
to differentiate between traces of airspeed” and glide slope deviation 
so another method was tried. 

The single channel 12 inch oscilloscope mentioned earlier was 
used in conjunction with a time sharing and blanking program. Thus 
a one amplifier scope was used to give the appearance of two independent 
and recognizable traces. The blanking circuit eliminated unwanted 
portions of the trace. A dot was chosen to represent airspeed deviation 
from trim where a "fast*' is above the datum and a "slow" is below. 



Glide slope information was displayed as a horizontal line. The 
line was beneath datum for a below glide slope signal and above 
for a "high-" This was opposite to the standard Instrument 
Landing System, "fly to the needle," display. Full scale on 
airspeed was four knots while full scale on glide slope deviation 
was eight feet. Photographs of the actual display are in Figure 
37 . 

The tasks required of the pilot were essentially the same 
as for the automatic systems. He was required to maintain glide 
slope and airspeed in the influence of gust disturbances. The 
additional task of responding to a "fly up 10 feet" command was 
also required. 

The pilot was thoroughly briefed before each task and 
allowed to practice until he was satisfied with his performance. 
He was told when and what kind of disturbance to expect on each 
test. The best of his efforts was retained. Thus the pilot 
was given every favorable chance to compete with the automatic 
system. 

The results of this Analog Simulation Program are presented 
in Chapter V. 
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CHAPTER V 



RESULTS AND DISCUSSION 



General 

Data from the analog simulation program in Chapter IV were 
obtained from the computer in the form of time histories. The 
traces were produced by a six channel Brush Recorder and are 
presented as Figures 38 through 72. Scales for the perturbations 
were kept uniform wherever possible for ease in comparison with 
the basic airplane. The basic airplane plus all automatic systems 
were calculated in time scale on the computer in order to speed 
up the solution. The traces, Figures 38 through 60, were run at 
10 mm per second, but in real time this is 1 mm per second. The 
manual simulation was accomplished in real time, and the records. 
Figures 61 through 72, were run at 5 ram per second in order to 
smooth out the appearance of the control inputs. 

Due to equipment malfunctions, the recording channels 
available were reduced to four by the time manual simulation was 
begun. Unfortunately the channels available were not grouped 
together, thus making it necessary to cut out each trace and fix 
the collection together with rubber cement. This was at best a 
tedious task. Nevertheless, a formidable amount of data was 
produced, and some manner of summarizing the important features 
was needed. 



46 



Measures of Effectiveness 



It was decided to group the data into responses to gusts and 
commands. The gusts were 1) 5 kt tail gust, 2) 5 kt up gust, 
and 3) combination of 1) and 2). The various commands were: 

1) step stabilator, 2) step spoiler, 3) one degree pitch, and 
4) fly up 10 feet. Thus grouped, there were obvious desirable 
quantities apparent. Such summaries of data are presented in 
Tables IV, V, and VI. 

The measures of effectiveness of a system in response to a 
gust disturbance were taken to be: 1) first time u = 0, 2) per- 
centage overshoot in u, 3) maximum h, and 4) value of h at 5 
seconds after disturbance. "First time u = 0" is a measure of 
the rise time of the system. This time is to be minimized but 
not at the expense of "Percentage overshoot in u."* The over- 
shoot is a measure of system damping. An arbitrary range of 
values from 15 to 25 per cent is considered optimum for over- 
shoot. It is evident that maximum glide slope deviation should 
be minimized. Of vital importance to the pilot is deviation 
from glide slope 5 seconds after disturbance. It takes about 
5 seconds to traverse the last 900 feet of glide path at a 
closure rate of 105 knots. The disturbance is assumed to occur 
at this point (Recall Figure IS) and could well mean the difference 
between a successful arrestment or a catastrophic collision with 
the flight deck ramp of the carrier. 

An evaluation of the responses to commands indicated the 
following measures of effectiveness: 1) time for h = 10 feet, 
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2) percentage overshoot in h, 3) maximum deviation in airspeed, 
and 4) peak values of(X • 

The maximum allowable deviation from glide slope at 900 feet 
out is 10 feet. This is the lower limit of the mirror landing 
system cone at this range. Any deviation lower than 10 feet 
will probably result in a ramp strike while 10 feet high could 
cause the airplane to "bolt," i.e., miss all the arresting wires. 

The time to eliminate this error obviously must be within the 
time envelope of 5 seconds mentioned above. Hence "time for h = 

10 feet" must be minimized. Again this time must not be minimized 
at the expense of overshoot since correction for a "high" could 
result in a "scooping out" at the ramp. The next measure was 
maximum deviation in airspeed; this is to be minimized. Excursions 
in angle of attack are also to be minimized. 

Gust Response 

Figures 38, 39, and 40 show the basic airplane response to 
gust disturbances. The phugoid oscillation is a readily apparent 
long period variation in pitch and airspeed while remains 

practically constant. The motion is lightly damped, and the 
period is 36 seconds. This result confirms the. well known assumption 
employed earlier. The short period is characterized by a heavily 
damped oscillation in(X . Here the period is six seconds and the 
motion is damped out after one cycle. The basic airplane results 
are important not only because they are used to compare the effects 
of adding control, but they also provide a means to check the 
analog model. 
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Reference 15 shows essentially the same results as shown 
in Figures 38 and 39, This favorable comparison was made even 
though Chance-Vought includes some non-linear effects in their 
analog model of the F-8. 

APC lightens the pilot’s load during a carrier landing 
approach by allowing him to concentrate more on glide slope 
and line up. The result is a more precise, and thus safer, carrier 
approach. APC is not and was not intended to be glide slope 
control. 

The performance of the APC was best determined from its 
gust response. Measures of effectiveness for the APC were taken 
from Reference 9 and are listed in Table IV. Here the APC shows 
the expected improvement over the basic airplane (configuration B) 
and the conventionally controlled airplane. Of interest is the 
case where the elevator is manually controlled but APC is used. 

Here it seems the pilot’s performance would have improved if he 
had not used the stick at all because the APC alone outperformed 
him in all measures. However, the APC plus pilot is better in 
performance than manual power compensation. 

Tail gusts cause the airplane to sink below glide slope while 
up gusts tend to make the airplane go high. Thus one test for the 
glide slope controller was gust performance. Summaries of glide 
slope controller effectiveness in gust conditions are presented 
in Tables V and VI. 

Inspection of the manual control records for the gust conditions 
yielded an average time lag for the pilot. The disturbance was 
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presented to the pilot by warning him, then suddenly illuminating 
his display scope. The time from disturbance to first control move- 
ment can be measured on the manual control record. The pilot's 
average time lag was found to be 0.5 seconds. 

Another test for the glide slope controllers was the response 
to a fly up command. This task was meant to test the system's 
recovery from a 10 foot low condition. These results are summarized 
in Tables IV and V. 

On the basis of the measures of effectiveness the best glide 
slope controllers were judged to be the DLC controllers. The DLC 
controllers were then compared with each other in Table V. The 
results showed that automatic DLC was superior to manual DLC, but 
the latter was a feasible means of control. 

The best DLC system was determined from gust and command 
performance. The best controller was, as expected, the DLC and 
APC system. 

The step commands. Figures 41 and 42, gave interesting though 
intuitive results. While step elevator initially caused the airplane 
to climb, the steady state result was a rate of descent; see Table V. 
The opposite was true for the suddenly retracted spoiler. A.n initial 
climb was followed by steady state rate of climb. Figure 41, Table V. 

Automatic glide slope control cannot be achieved by means of 
EGSC without the use of an APC. In fact, an}^ control of glide slope, 
automatic or not, is impossible with elevator control only. The stick 
and throttle must be skillfully coordinated in order to maintain or 
recapture glide slope. Thus the automatic systems currently being 
tested by the Nav’y (7) are ineffective without a functioning APC. 
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Figures 42 and 60 substantiated this intuition. Step stabilator 
for the intention of climbing results in steady state descent 
due to the increase in drag. A Pitch Control System also yielded 
a rate of descent after nose up pitch was commanded. Hence the 
glide slope controller using pitch control without APC initially 
corrected in the proper sense for a low, but the steady state 
response was improper. As the nose came up, airspeed decreased 
due to increased drag and decreased kinetic energy. The steady 
state value ofQ^became greater than pitch angle, 0. Hence from 
Equation 9, a negative glide slope angle resulted. 

Figure 60 also emphasizes the importance of the integration 
term used in all the other controllers in this context. Recall 
that the integration term was used to wash out steady state error. 
The PCS, Figure 18, did not incorporate this feature} hence, a steady 
state error in © appears in Figure 60. One degree of pitch angle 
was commanded, but only 0.6 degree was produced. 
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CHAPTER VI 



CONCLUSIONS AND RECOMMENDATIONS 



Cone lusions 

This study encompassed three areas of investigation: 

1) implementation of wing lift spoilers as a DLC system, 

2) automatic control of glide slope, and 3) the analog 
simulation program. 

Conclusions from the DLC implementation are as follows. 

First, wind tunnel data, preferably full scale, are needed 
when dealing with wing lift spoilers. Second, spoilers provide 
a means of rapidly increasing lift while decreasing drag. 

In the study of automatic glide slope control the following 
conclusions were made: First, not only position feedback but 

rate feedback as well, must be used in automatic glide slope 
control. Second, increasing the gain of the outer loop of the 
automatic glide slope controllers results in decreasing the 
system damping. Third, recourse to an analog computer for 
determination of gain constants is easier than the root locus 
method. 

First, the analog simulation program indicated that a DLC 
system using spoilers could be used for glide slope control 
without an APC . Second, an elevator glide slope coupler (EGSC) 
system incorporating an APC was not as effective as a DLC or 
APC system. Third, EGSC system cannot control glide slope 
automatically without a functioning APC. Fourth, a man controls 
glide slope remarkably well but is not as effective as the automatic 
sys terns . 
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Recommendations 



The work involved, during this investigation, in gathering 
data on spoiler effectiveness was almost insurmountable. Data 
were at best sparse and inconsistent. It is therefore recommended 
that a meaningful investigation be made of the characteristics of 
spoiler systems. 

The hybrid computer which is available at the Naval Post- 
graduate School opens the door for a wealth of simulation studies. 
The high speed digital retrieval of data will allow simulation to 
be accomplished using non-linear equations of motion and random 
disturbances . 
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TABLE I 



DIMENSIONAL STABILITY BERIVATI^TIS 
FOR THE F-8 IN APPROACH CONFIGURATION 



DATA: 



S = Wing area 
w 

= Thrust arm 
= Tail arm 
c = Mean chord 
B = Moment of inertia 
U = 1.3 X stall speed 
eg = Center of gravity 
W = Gross weight 
(Q = Pitch reference 
= Thrust angle 



= 375 ft^ 

= - .437 ft 
= 14.08 ft 
= 11.8 ft 
= 96000 slug ft^ 
= 234 ft/sec 
= ,24c 
= 22000 lb 
=8,1 deg 
= .85 deg 



SYMBOL' 


DEFINITION 


DIMENSIONALIZER 


VALUE 


X =: 

U 


1 C>X _ 
C)UL ■ 




= - o 060 1/sec 


X 

w 


-L ^ - 
rn dto 




= -.01419 1/sec 




_L SX - 
rn ■ 




= 2.48 ft/sec^ 




X dX - 




2 

= - .1.64 ft/sec rad 


It 

H 

X 


1 ax . 

§T - 




= .00145 ft/lb sec^ 


Z r 

u 


X .aZ . 

rn 


^ (-X 


= - .2655 1/sec 
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TABLE I (CONTINUED) 



SYMBOL DEFINITION DIMENSIONALIZER VALUE 



z = 

w 


m 


Z)V 




= - .4265 


1/sec 




ZE = 


1 

m 


as 


= 


= - 85 ft/sec^ 




= 


1 

m 




■ ®‘(-q 


= - 19.1 


f t /sec*" 


} 

“rad 


S = 


1 

m 


dT 


m 


= - 2.17 


X lO'^ 


f t/sec^rad 


M = 
u 


i_ 

& 


m 

aoL 


- pgy(-Gz,'^ 


= .000185 


1/sec 


ft 


M = 
w 


1 

b 


auo 


II 

o 


= - .004858 l/s( 


ic ft 


M. = 
w 


1 

B 


3H 

Sli> 




= - 1.772 


-1 

X 10 


• i/ft 


M = 

q 


1 


ad 

H 


■ ffS'fc",) 


= - .3384 


1/sec 


rad 




1 

& 


ad 




= - 2.25 


1/rad j 


2 

56C 


“x ' 


i 

5 


M 

ai 


= ?2Lx 

B 


= - 4.55 


X lo”^ 


1/lb sec 
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TABLE II 



ANALOG POT SETTINGS 



Pot 


Represents 


Setting 


Ga in 


1 


• 060 s-vX / 3 .U, 


.1011 


1 


2 


, 332 SiX / 


,0116 


1 


3 


31.85 au/a^ 


.111 


1 


4 


2.5 ao./as 


.0051 


1 


5 


.00145 au/a-t 


.0292 


10 


6 




.593 


1 


7 


,00113 aj^ /a^^ 


.547 


1 


8 


.4265 a^/a^ 


.4265 


1 


9 


.01935 aj(/ae 


.01935 


1 


10 


a^/ae 


.1 


10 


11 


. 08 15 aj( / a^i^ 


.00815 


10 


12 


,353 a^/a^ 


.202 


1 


13 


a© /a^ 


.286 


10 


14 


1.85 X 10 ^ a^/a^ 


.0312 


1 


15 


.338 a^/a© 


. 118 


1 


16 


.0415 a^/a^ 


.0415 


1 


17 


1 . 14 a^ /a<;t 


.0398 


10 


18 


4.55 X 10 ^ a^/aT 


.0091 


10 


19 


2.25 


.0786 


10 


20 




Various 


1 


21 




Various 


1 


22 


U ay^/a^ 


.409 


1 
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23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 



TABLE II (CONTINUED) 



Represents 


Setting 


Gain 


U Kj:, aj, /a, 


.300 


10 


ag/ay, 


.500 


10 




.001 


1 


he 


Various 


1 


Ku, a-r/a^ 


.337 


10 


a-f / a^_ 


.0336 


1 


a-T / a^^ 


.1745 


1 


Te 


.8711 


1 


Te 


.8711 


1 


U a^/a, 


.409 


1 


U Kv,a^/a^ 


.250 


1 


he 


Various 


1 


K'n a,^/ah 


.300 


1 


Ky^Kj_a^ /ay^ 


.0025 


1 


u 

CD 


Various 


1 


Kea^ / 


.13 


10 


Ke ae/a^ 


.677 


1 
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TABLE II 



Parameter 

u 

U , Uc , Ucj 

cx,o^,<(.e, 

Q 

T 

h, he 

S 



ANALOG SCALE FACTORS 



Maximum Value 
50 ft/sec^ 



Scale Factor (a^) 
2 

2 volts/ft/sec 



84.3 ft/sec 
5 kts 



1.189 volts/ft/sec 
2 volts /kt 



10 deg 
.1745 rad 



10 volts/deg 
573 volts/rad 



10 deg/sec 
.1745 rad/sec 



10 volts/deg/sec 
573 volts/rad/sec 



.5 rad/sec^ 



200 volts /rad/sec" 



1000 lbs 



.01 volts/lb 



100 ft 



1 volt/ft 



.Ic 



1000 volts/c 
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TABLE IV 



EFFECTIVENESS OF APC COMPARED 
WITH BASIC AIRPLANE AND MANUAL CONTROL 



Dis turbance 
Figure Number 
Thrust Control Type 
Glide Slope Control System 
Glide Slope Control Type 
First Time u = 0, Sec. 

Time for -l<u<^l, Sec. 
u Overshoot, % 

Maximum Thrust Change, Lbs, 



Disturbance 
Figure Number 
Thrust Control Type 
Glide Slope Control System 
Glide Slope Control Type 
First Time u = 0, Sec. 

Time for -Ku<,l, Sec. 

u Overshoot, 7o 

Maximum Thrust Change, Lbs. 



Disturbance 
Figure Number 
Thrust Control Type 
Glide Slope Control System 
Glide Slope Control Type 
First Time u = 0, Sec. 

Time for -l<u<l. Sec. 
Maximum Thrust Change, Lbs. 



5 Knot 


Tail 


Gust 




38 


47 


63 


66 


B 


AUTO 


MAN 


AUTO 


B 


B 


EGSC 


EGSC 


B 


B 


MAN 


MAN 


8 


2.5 


6.5 


3 


102 


2 


5.2 


2.5 


85 


33 


50 


40 


— 


1800 


-600 


2600 


5 Knot 


Up Gust 




39 


48 


62 


65 


B 


AUTO 


MAN 


AUTO 


B 


B 


EGSC 


EGSC 


B 


B 


MAN 


MAN 


12 


20 


17 


8 


10 


0 


11 


5 


80 


0 


0 


0 


— 


300 


-600 


-650 


Combination 


Gust 




40 


49 






B 


AUTO 






B 


B 






B 


B 






9 


3 






102 


2.5 






— 


2000 
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TABLE V 



SUMMARY OF GLIDE SLOPE 
CONTROL EFFECTIVENESS 



Dis turbance 


5 Knot 


Tail 


Gust 






Figure Number 


38 


43 


50 


55 


63 


Thrust Control Type 


B 


B 


AUTO 


AUTO 


K4N 


Glide Slope Control System 


B 


DLC 


DLC 


EGSC 


EGSC 


Glide Slope Control Type 


B 


AUTO 


AUTO 


AUTO 


MA.N 


First Time 0=0, Sec. 


8 


60 


3 


3 


6.5 


u Overshoot , % 


85 


0 


20 


14 


50 


h Maximum, Ft. 


-115 


-6 


-4 


_7 


3 


h at 5 Sec . , Ft , 


-20 


-6 


-4 


-6 


-2 


Disturbance 


5 Knot 


Up Gust 






Figure Number 


39 


44 


51 


56 


62 


Thrust Control Type 


B 


B 


AUTO 


AUTO 


MAN 


Glide Slope Control System 


B 


DLC 


DLC 


EGSC 


EGSC 


Glide Slope Control Type 


B 


AUTO 


AUTO 


AUTO 


MAN 


First Time u = 0, Sec, 


14 


7 


5 


15 


17 


u Overshoot , % 


80 


50 


— 


0 


0 


h Maximum, Ft. 


27 


3 


3 


5 


5.5 


h at 5 Sec . , Ft . 


15 


0 


4 


5 


1 


Disturbance 


Combination 


Gust 






Figure Number 


48 


45 


52 


57 




Thrust Control Type 


B 


B 


AUTO 


AUTO 




Glide Slope Control System 


B 


DLC 


DLC 


EGSC 




Glide Slope Control Type 


B 


AUTO 


AUTO 


AUTO 




First Time u = 0, Sec. 


9 


60 


3 


3 




u Overshoot, % 


80 


0 


20 


20 




h Maximum, Ft. 


-105 


-6 


-4 


-3 




h at 5 Sec . , Ft . 


-5 


5 


-3 


-1 
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TABLE V (CONTINUED) 










— — — 


— 




Command 


Fly Up 


10 Feet 








Figure Number 


46 


53 


58 


61 


64 


Thrust Control Type 


B 


AUTO 


AUTO 


MAN 


AUTO 


Glide Slope Control System 


DLC 


DLC 


EGSC 


EGSC 


EGSC 


Glide Slope Control Type 


AUTO 


AUTO 


AUTO 


MAN 


MAN 


Time to 10 Ft., Sec. 


3.5 


3.8 


7 


5.3 


5 


h Overshoot, % 


20 


25 


25 


25 


20 


u Maximum, Kts . 


— 


.75 


-.8 


-1 


-.6 


(XFeak, Deg, 


1.1 


.5 


.5 


— 


— 


Command 


Low/Slow (10 Ft 


. /5 Knots) 




Figure Number 


54 


59 








Thrust Control Type 


AUTO 


AUTO 








Glide Slope Control System 


DLC 


EGSC 








Glide Slope Control Type 


AUTO 


AUTO 








Time to 10 Ft., Sec. 


6 


12 








h Overshoot, % 


25 


30 








u Maximum, Kts. 


.3 


-1 








CXPeak, Deg. 


-.5 


1 








Command 


Step Control Input 






Figure Number 


42 


41 


60 






Thrust Control Type 


B 


B 


B 






Glide Slope Control System 


EGSC 


DLC 


PCS 






Type Input 


1 


s 


G 






^ Initial 


+ 


+ 


+ 






^Steady, Deg. 


-1.75 


2 


-2 






©Peak, Deg. 


2.4 


1.3 


.9 






(XPeak, Deg. 


1.2 


-.2 


2.1 
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TABLE VI 



COMPARISON OF EFFECTIVENESS OF 
MANUAL AND AUTOMATIC DLC SYSTEMS 



Disturbance 
Figure Number 
Thrust Control Type 
Glide Slope Control System 
Glide Slope Control Type 
First Time u = 0, Sec. 
u Overshoot, % 
h Maximum, Ft. 
h at 5 Sec . , Ft . 



Disturbance 
Figure Number 
Thrust Control Type 
Glide Slope Control System 
Glide Slope Control Type 
First Time u = 0, Seo. 
u Overshoot, % 
h Maximum, Ft. 
h at 5 Sec . , Ft . 



Command 
Figure Number 
Thrust Control Type 
Glide Slope Control System 
Glide Slope Control Type 
Time to 10 Ft . , Sec . 
h Overshoot, % 
u Maximum, Kts. 

0(Peak, Deg. 



5 Knot 


Tail Gust 




69 


72 


43 


50 


MAN 


AUTO 


3 


AUTO 


DLC 


DLC 


DLC 


DLC 


MAN 


MAN 


AUTO 


AUTO 


1 :i 


3 


60 


3 


50 


20 


0 


20 


7.5 


-7 


-6 


-4 


-5 


-7 


-6 


-4 


5 Knot 


Up Gust 






68 


71 


44 


51 


MAN 


AUTO 


B 


AUTO 


ELC 


DLC 


DLC 


DLC 


MAN 


MAN 


AUTO 


AUTO 


13 


6 


7 


5 


25 


0 


50 


0 


5 


4 


3 


3 


-1 


0 


0 


4 


Fly Up 


10 Feet 






67 


70 


46 


53 


MAN 


AUTO 


B 


AUTO 


DLC 


DLC 


DLC 


DLC 


MAN 


MAN 


AUTO 


AUTO 


11 


6 


3.5 


3.8 


15 


30 


20 


25 


-2 


■1 


-1 


.75 


— 


— 


±•5 


+.5 


__ , 


■ - 
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TABLE VII 



COMPONENT TRANSFER FUNCTIONS 



I 1 Ih 

I 1 ’1,154 lb 


- _i ■?. JC£)A. 
6a ' rad 


a - -'o.ooo 


1®-- -0.667 sec 


_a. -4.5500-'“ lai 

T - asO. 8064 0,281 lb 


I- i 


Os 4 J ft 


G . -(2.25i+a867Wa<l 

y ,4(3*'’0.60b.4’l.Z8l) rod 


U . 0,001 45 ,a ft 

T ' 4N0.06^ + 0.036 seclb 


e . -(2.25/ltO.&67Uad 
Y ^^-^0. 806 -r Z8 1 ro4 


-0.0073 ^ 


(X _ - (6.3G^\o.lZZ^)ra4 
8 - .atAQ&Qbdtl 26\) 


- 0 . 003(1 ^ 

he \ A J {i 


Q . 0.01425+0.400 rad 

S ’i(s‘4),80b^428l) 


h - 1 Sec 

h 


2 .O.56/.0.1 36^.0,409 red 



65 



Fuselage 




FIGURE 1 

F8 WING PLANFORM 

b =35 6 FT 

X= .28 



S = 375 FT 
w 



AR = 3.3 



6b 



Fuselage 




FIGURE 2 

NACA RM A54H26 MODEL 2 
S = 312.5 FT.* b = 30.6 FT 



.7c 



AR =3.0 X = .4 



6 / 




VARIOUS SPOILER DEFLECTIONS 
NACA RM A54H26 MODEL 2 
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A Cl vs. S for various a 

NACA RM A54H26 MODEL 2 



69 



1 . 2 - 



1 . 0 - 



.8 ■ 



.6 • 



APPROACH C[_ =0.905 



. 2 ^ 




0 



1 

I 



4 8 12 16 20 

a DEG 

w 

FIGURE 5 

F8 - Cl vs. 



EFFECT OF SPOILER DEFLECTION (EST.) 



IS 



2.4 



I 
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FIGURE 6 

Cp vs. a 

VARIOUS SPOILER DEFLECTIONS 
NACA RM A54H26 MODEL 2 
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1 1 1 

12 13 14 15 16 



a, DEG. 

FIGURE 7 

C (3 vs. a IN DLC RANGE 

VARIOUS SPOILER DEFLECTIONS 
NACA RM A54H26 MODEL 2 
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□ a . 14* 

o 15* 



. 0 2 - 



.0 1 - 




FIGURE 8 

aC, vs. S for various a 
NACA RM A54H26 MODEL 2 
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. 3 5- 
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I 5- 
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RANGE 



.. . , 

0 4 8 12 16 20 

a, DEG. 

FIGURE 9 

F 8 C vs. a 
D w 

EFFECT OF SPOILER DEFLECTIONS (EST.) 
C = -.07 
DS 
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approx. F8-CG (?.27 MGC 




FIGURE 10 

^^cg' 

EFFECT OF SPOILER DEFLECTION 
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FIGURE 11 
AXIS SYSTEM 

ARROWS SHOW POSITIVE DIRECTION 
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BLOCK DIAGRAM - ARC 
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FIGURE 13 BLOCK DIAGRAM -DLC (AUTO) 
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FIGURE 14 BLOCK DIAGRAM - EGSC (AUTO) 
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FIGURE 16 

ROOT LOCUS PLOT 
DLC FEEDBACK LOOP 
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REAL 
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DLC OUTER LOOP 
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FIGURE 19 

APPROACH GEOMETRY 



I 
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FIGURE 20 ANALOG DIAGRAM - BASIC AIRFRAME 
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ANALOG DIAGRAM — DLC (AUTO) 
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FIGURE 22 




fO 

CM 

lU 

q: 

3 

S2 

u. 



87 



ANALOG DIAGRAM — EG SC (AUTO) 
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ANALOG DIAGRAM - PITCH CONTROL SYSTEM 



FIGURE 25 

ARC GAIN OPTIMIZATION 
U RESPONSE TO 5 KT TAIL GUST 
Ku= 400 KuKI* VARIOUS Ka = 10,000 
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FIGURE 26 

ARC GAIN OPTIMIZATION 
U RESPONSE TO 5KT TAIL GUST 
Ku = VARIOUS KuKi-,12 Ka = 10,000 
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FIGURE 27 

ARC GAIN OPTIMIZATION 
U RESPONSE TO 5KT TAIL GUST 
Ku =400 Ku Kj = .12 Kq = Various 
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Id - M 
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FIGURE 28 

DLC GAIN OPTIMIZATION 
RESPONSE TO 10 FT. COMMAND 
^r0073 Kh = VARIOUS KhKj = IO' 




Id - ^ 
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RESPONSE TO 10 FT. COMMAND 
= VARIOUS Kh = .005 Kh K| = 10’® 




I I I 



Id - ^ 
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DLC GAIN OPTIMIZATION 
RESPONSE TO 10 FT. COMMAND 
= -.0073 Kh = .005 KhKi = VARIOUS 
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6 r 



Kh=8.66 X lO'"^ rad /ft 
Kh« 5.2 X 10"^ 




FIGURE 33 

ARC + EGSC GAIN OPTIMIZATION 

h RESPONSE TO 10 FT COMMAND 

Kh = VARIOUS Kfi = 4.4 X lO""^ Kh K| = I.73 X lO"® 
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GENERAL ARRANGEMENT 




FIGURE 35 

THROTTLE AND CONTROL STICK 
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FIGURE 37a 

PILOT'S VISUAL DISPLAY SHOWING 
ON SPEED 2 FEET BELOW GLIDE SLOPE 




FIGURE 37b 

PILOT'S VISUAL DISPLAY SHOWING 
2 KNOTS SLOW ON GLIDE SLOPE 
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FIGURE 4l 
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“»«'«»• •• FIGURE 47 
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u -Krs 



-5 



—5 SEC 






*4000 

T -lbs 



RESPONSE TO 5 KT UP GUST 
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T-'LBS 



FIGURE i^.9 
AIRPLANE + APG 

RESPONSE TO 5 KT TAIL + 5 KT UP GUS 




‘4000 
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*4000 
T -LBS 



AIRPLANE + DLC + APC 
RESPONSE TO 5 KT UP GUST 



-4000 
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5 SEC 
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AIRPLANE + DLC + APC 
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FIGURE 53 

airplane + JD(LC + APC 
RESPONSE TO 10 FT h COMMAND 
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FIGURE 56 

AIRPLANE + EGSC + APG 
RESPONSE TO 5 KT UP GUST 
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FIGURE St 

AIRPLANE + EGSC + APC 
RESPONSE TO 5 KT TAIL GUST 
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FIGURE 58 

AIRPLAME + EGSC + APC 
RESPONSE TO 10 FT h COMMAND 
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FIGURE 59 

AIRPLANE + EGSC + APC 
RESPONSE TO 5 KT TAIL GUST 
AND 10 FT h COMMAND 
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Q-D£G 



-1 






FIGURE 60 



AIRPLANE + PITCH CONTROL SYS'TEM 
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FIGURE 61 

MANUAL THROTTLE + MANUAL STABILATOR 
RESPONSE TO 10 FT h COMMAND 
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FIGURE 62 

MANUAL THROTTLE + MANUAL STABILATOR 
RESPONSE TO 5 KT UP GUST 



126 




FIGURE 63 

MANUAL THROTTLE + MANUAL STABILATOR 
RESPONSE TO ^ KT TAIL GUST 
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FIGURE 64 

APG + MANUAL STABILATOR 
RESPONSE TO 10 FT h C01H>1AND 
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FIGURE 65 

APC + l^NUAL STABILATOR 
RESPONSE TO 5 KT UP GUST 





FIGURE 66 

APC + MARUAL STABILATOR 
RESPONSE TO 5 KT TAIL GUST 
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FIGURE 6/ 

MANUAL THROTTLE + MANUAL SPOILERS 
RESPONSE TO 10 FT h COMMAND 
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FIGURE 68 

MANUAL THROTTLE + I4ANUAL SPOILERS 
RESPONSE TO 5 KT UP GUST 
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PIG'JH3 69 

iiAUU.lL TinO^’TLE + MAiniAL GPOILURS 
R9:3PCir3E TO 5 KT TAIL Gj3T 
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FIGURE 70 

APG + MANUAL SPOILERS 
RESPONSE TO 10 FT h COMMAND 
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FIGURE 71 

AFC ■* MANUAL SPOILERS 
RESPONSE TO 5 KT UP GUST 
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FIGURE 72 

APC + MANUAL SPOILERS 
RESPONSE TO 5 KT TAIL GUST 



136 



APPENDIX 



DETERMINATION OF LOOP GAINS FOR AUTOMATIC SYSTEMS 



APC 

For the purpose of determining approximate loop gains, the 
inner loop of the APC system was first analyzed independently 
of the rest of the system. Referring to Figure 12, a simplified 
inner loop was constructed. Figure 73. 
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FIGURE 73 

Inner Loop, APC System 

From Figure 73, using the short period assumption for the QC 



loop. 



T. = 



T 



A 

TB=Tc 

TZ T1-+ lju. = Tc 



Combining these three equations yielded 



Ta tlTu -fTB 



I 

A 



Replacing A and B by their respective transfer functions from 
Table VIII and neglecting the engine time lag in order to reduce 
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the characteristic equation to an easily handled second order 



function, 

T _ A^O-SObi -> I.2 &1 

1,281.4.55^0-“ K. 

The characteristic equation is then, 

2)40.80!o/i 1 1.2&1 ^ 4.55 - 0 

from which, 

Z 0.606 

and 

wJ. = l.28U4.55»lO'‘'Ko< 

For a damping ratio of ^ = 0,5, simultaneous solution of these 
two equations yielded, 

Kc< ~ - 1. 485 AO i bs/ro.di 

The inner loop of the APC system was then reduced to the single 
transfer function below. 

T . / j% Q.&OG ^ ^ 1.Z81 

/S' t0.80G i) t0.05 

For the analysis of the overall loop, ~|"^j^was neglected since 
it is small compared to . A simplified loop, derived from 
Figure 12, was used for the analysis. The simplified loop is shown 
in Figure 74. 

, . i I u 

^ Uv l.'j. I 

FIGURE 74 

Outer Loop, APC System 
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Using the phugoid assumption for , the simplified transfer 
function for- the APC was derived in the following manner, 



a 



_ LL 



e - 



c 



Lie- U -■ Lie 

The result of combining these two equations was. 



LI - 1 

Uo' 1 4- C 

Insertion of values for the transfer functions in C from Table 
VIII and from the inner loop analysis gave 

o. .QQ145Ka A-^.QQli7K^ ^ .QO\&eA 

aJ AM8-G6..00l45Kai^\(.745h00ll7Ka)4M,0G&t,0O186Ka)i+.0Z3 



from vhich the characteristic equation is 

^M.8ia6,.OOI45U4M743-.OOil7K.»)AM.068-,OOl8fc Ui ,.0Z3 -0 

A digital computer program was used to make a root locus plot 
of the characteristic equation, Figure 15. From this plot, for 
a damping ratio of 0.5, an estimate for sec/ft was 

found . 



DLC 

The DLC system was analyzed in a manner similar to the above. 
The short period assumption was used throughout and the integration 
loop was neglected. The simplified inner loop block diagram, taken 
from Figure 13, is shown in Figure 75. 
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FIGURE 75 

Inner Loop, DLC System 



From Figure 75, the following analysis was made 

S ^ 

Sh= Pit 
Eh-v Sh= £} 

These three equations were combined to form 

Jq_ - D 
Sk' 1 -DE 

Incorporation of numerical values from Table VIII for D and E 
resulted in, 

Vi - ^^4(. .154 ^ -t. 41 

S ■ •t(,806-84.5Kt,)A’ 4(1.1811 3l.4KhU -96K^ 

From this, the inner loop characteristic is 

AM.80fe-84.3KnlA’4(l.7,8l-5l.4Ki,)a - %K<>=0 

The root locus plot for this characteristic equation is found in 
Figure 16. For a damping ratio of 0.8, the resulting was 
found to be -0.003 sec/ft. The resulting transfer function for 
the inner loop is 

W . Zm.'5bA^^ 4 . 4 '') 

Sh ^ , Z85 

The simplified block diagram for the DLC system was derived 
from Figure 13 and is shown in Figure 76. 
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FIGURE 76 

Outer Loop, DLC System 



Analysis of Figure 76 yielded 



from which 



he= 

hc-h= he 

F 



h - 

he 



t r 



The transfer function for the DLC system is then, 

2 34Kh(.36^%.\!)4z>-t.4n 






/.Sl,0b,!l^58+a4,3K„')4M.268t51.4K»U +96K,^ 



from which the characteristic equation is 

2)M.06A^(l.!)8t84.5KJ4^(.28a + 51.4Kh)At%Kh -0 

Root locus analysis of this function yielded a value of 0.005 for 

Ko . 



PCS 

Figure 77 shows the block diagram for the simple Pitch Control 
System. This system was used to demonstrate the lack of feasibility 
of conventional elevator, without speed control, as a means of 
automatic glide slope control. 
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FIGURE 77 

Simplified PCS System 



Using the short period assumption for , the following 
analysis was made 









'1 





‘o( ■ 


1 \j 


11 
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.G . 


LMi. 


! { J 



from which 



and 



- - (2.2 5 b 

)/ ■ 2(2)S 8062 t 1.261) 

Q - 23 ^ . -( 2.2 5 A + , 867) 

'q ■ " 2^)^t,806Z)^ 1.Z81 

Referring to Figure 77, the inner loop transfer function is 

A - H 



I 3-HL 
-(2.25/>f.8GT^ 



from which 

Q _ 

■%" 4*^(.30to42.25foU ,(1.28i-.Su 

The characteristic equation for the inner loop is 

L8Ck.-> t 2 •/ jKo) .1) - ( l./ibl' Gc K.aG '■ G 
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From the characteristic equation, 

-- l.28\ - .867Ke 
806-2. 25 Ks 

When these two equations were combined, the result was a quadratic 
equation in |<^q . 



Solution of this equation yielded, 

K© = 989 j , 667 ( <i>rf\ef\sioniess^ 

Returning to the characteristic equation, it was seen that, for 
the system to be stable, 

.806-2.25Ki>0 

and 

U8I -,867Ke>0 



These two inequalities imply that 

Ko< ,36 

Therefore the correct value for is -0.677* The resulting 
transfer function for the (5 loop is, 





-( 2.75 2 ^+. 867 ) 






G . 


- 12 .25^4 .867) 






The simplified block diagram for the PCS outer loop is shown 
in Figure 78. 
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FIGURE 78 

Simplified Block Diagram, PCS Outer Loop 
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Analysis of this loop yielded the result, 

0 

0c' > 

from which 

^ - -(2.25 Ke J i . 8 67 Kb 'i 

Qc 2)^ + 2, 3264Sn, 858-2. 25Ko'l/i-. 867 Ke 

The characteristic equation is 

2.5ZG 1.558-2.25 K©V<^ -.807 Ke=-0 

The root locus analysis of this characteristic function resulted 
in the choice of )<^^= -1.3 for a damping ratio of = 0.5. The 
root locus plot for the outer loop of the PCS system is shown in 
Figure 79. 
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